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Abstract. Often the current mode coupling theory(MCT) of glass transitions 
is compared with mean field theories. We explore this possible correspondence. 
After showing a simple-minded derivation of MCT with some difficulties we give 
a concise account of our toy model developed to gain more insight into MCT. We 
then reduce this toy model by adiabatically eliminating rapidly varying velocity- 
like variables to obtain a Fokker-Planck equation for the slowly varying density- 
like variables where diffusion matrix can be singular. This gives a room for non- 
ergodic stationary solutions of the above equation. 



1. Introduction 

Our understanding of phase transition starts from the famous thesis of van der Waals of 
1873 where complex effects of intermolecular interactions were put into two parameters 
often denoted as b and a representing, respectively, the repulsive and attractive parts 
of intermolecular forces. This theory was transcribed into magnetism resulting in the 
Weiss theory of ferromagnetism which contains a single parameter measuring strength 
of the molecular field. These mean field theories were remarkably successful and 
provided a good beginning of phase transition theory. Shortcomings of these theories 
were soon noticed, especially after Onsager's exact analysis of the two-dimensional 
Ising model. Thus the basis of the mean field theories were examined and exactly 
solvable model was constructed which yielded mean field theory results The 
characteristics of this model is that the attractive part of intermolecular potential 
has an infinite range after a certain limiting procedure, which suppress fluctuation 
effects responsible for deviations from the mean field behavior. Efforts to incorporate 
neglected fluctuation effects led to our current day understanding of phase transition, 
in particular, critical phenomena. 

In comparison, the case of structual glass is much less clear-cut Q]. First there 
are still uncertainties about nature of the glass transition. It is not clear whether there 
is a real transition or just a cross-over. Eventhough one assumes a genuine transition, 
opinions differ among those who think that there is an underlying thermodynamic 
transition and those who believe that the transition is merely a kinetic one. The case 
of spin glass is much more well understood, at least on the mean field level where we 
know the existence of a genuine thermodynamic phase transition 

Under such circumstances it is quite significant that the first principle theory of 
structural glass transition, that is, the MCT, was proposed and succeeded in explaining 
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some aspects of glass transitions y, Since the theory is still rather crude, it is fair 
to regard the status of this theory as similar to that of van der Waals and Weiss 
theories mentioned above. As in the case of these theories the current MCT is beset 
with serious difficulties. These basically come from the fact that the MCT formalism 
was originally developed for critical phenomena focussing on very large length scales 
reaching to thousands of Angstroms. Applications to glasses have difficulties due to 
its primarily short length scales and due to its own peculiarities. These difficulties are: 
(a) The factorization approximation which replaces the four-body time corrrelation 
functions by the product of two-body time correlation functions is essential to obtain 
the self-consistent MCT equation. This is especially uncontrolled at short length 
scales of at most 10-20 Angstroms. See the following section, (b) The idealized MCT 
predicts a sharp dynamic transition to a nonergodic state at a certain temperature. 
But MCT does not provide any information on the nature of this nonergodic state, 
(c) The physical picture of the so called hopping processes in an extended version of 
MCT is stiff lacking. 

It should then be an urgent task for further progress to clarify the bases of the 
MCT. Motivated by this desire we constructed a toy model having the following three 
features 

• reversible mode coupling mechanism 

• trivial statics 

• mean-field type so that the model can be exactly solvable. 

This toy model is distinguished from other toy models for glasses in that it closely 
mimics the MCT g |. 



2. MCT 



Here we present an over-simplified derivation of a self-consistent MCT equation of the 
density-density time correlation function obtained first in We start from the 

following hydrodynamics- like continuum equation Q : 

d 

m— p(r,t) = -V-j(r,i) (1) 

|j(r,t)=f(r,t) + ... (2) 

Here m is the mass of a fiuid molecule, p(r,t) is the number density, and j(r,t) is 
the momentum density. The ellipsis in (||) contains the terms second order in the 
momentum density, dissipative terms and thermal noise terms, which do not play a 
role here and will be dropped hereafter. 

Here we assume the existence of a free energy density functional H{{p}). Then 
the body force density f(r, i) is the number density times a force on a test particle 
of the same kind. The latter is negative of the gradient of the infinitesimal variation 
of the free energy density functional against infinitesimal density change, that is, 
—'VSH{{p})/6p{r). Therefore we find (Here time arguments are omitted.) 

f(r) = -pir)ySHi{p})/Sp{r) (3) 

Now, the exact form for H{{p}) is unknown and various approximate forms are 
proposed. For the purpose of deriving the MCT equation, it suffices to use the popular 
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Ramakrishnan-Yussouf form despite its shortcomings 
H{{p}) - fcsT j dvp{v 



ln(^)-l 



-\kBT JdrJ dr'c(|r-r'|)(p(r)-po)(p(r')-po) (4) 

Here po is the density of the reference uniform liquid and c(|r — r'|) is the direct 
correlation function |^ of the reference liquid. The Fourier transform c(k) of c(r) is 
connected with the static structure factor of reference liquid SCk) through 

c(k) = - ^(k)-i (5) 

The fact that the direct correlation function appears here is important because this is 
the only place in this theory where the short range correlation central to any liquid 
theory is incorporated. 

The next step is to split f(r,i) into terms linear and quadratic in the density 
difference Sp{r, t) = p{r, t) — pq as 

f(r,<) = f'(r,<)+r'(r,i) (6) 
f'(r,t) = -Vp(r,t) (7) 

r\r,t) = kBT j dT'c{\v-v'\)5p{T,t)V'5p{v',t) (8) 

where p(r, i) is the local pressure correct up to 5p{v,t). Combination of and 
sind (H) tells us that the f' produces only linear density oscillations with 
constant wave vectors and does not contribute to freezing. On the other hand, f"'(r, t) 
is the sum of numerous terms oscillating with different frequencies, which, on the 
whole, look quite irregular. This fact also makes it hopeless to try to find solutions to 
these equations. 

However, we are not interested in individual solutions but only some statistical 
properties of them, which are also measurable quantities. Among such quantities 
the most attention is paid to the density-density time correlation function, whose 
normalized form is defined in terms of /Ck(i), the Fourier transform of Sp{r,t), by 

, ... _ < Pkft)p-k(O) > 
'^'^'^ = < Pk(0)p-k(0) > 

In obtaining this quantity from (Q) and (||) we can regard f"'(r,t) as a kind of 
random force familiar in the Langevin equation of Brownian motion [p^ . Then the 
equation that determines (/)fe(t) requires knowledge of a memory kernel which is the 
time correlation function of f^' (t) , the Fourier transform of f (r, t) . 
Explicitly the equation for (f)k{t) turns out to be 

'^ = -^lMt)-jysM,it-s)'-^ (10) 



where ftk = k^/ksT / Sk is the frequency of the local density oscillation in liquids. 
Aik{t) is the memory kernel given by 

■^^(*) = r^™7(/k'W/-k(0)) (11) 

mpoKB-i k^V 

where = ik • fk"' and V the system volume. Since this correlation function 
involves products of four density fluctuations which are impossible to deal with 
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directly, this is factorized into products of two density-density corrrelation functions. 
This produces the feedback mechanism responsible for freezing. The resulting self- 
consistent equation for the density-density correlation function is the same as that 
given for the first time inQ. In the above simple derivation we have side-stepped the 
fact that the time dependence of f"'(r, t) in the memory kernel is in fact governed by 
the "projected" dynamics in the sense of Zwanzig-Mori formalism. The correct but 
more involved derivation starts from a Fokker-Planck type equation . Alternative 
simpler derivation explicitly relies on the strong assumption that density fluctuations 
at various times obey Gaussian statistics . 

The transition to non-ergodic states in this theory is driven by the nonlinear force 
term (^), which, in turn arises from the quadratic (or harmonic) term of (^). The first 
term there containing a logarithm is just for ideal gas. Therefore, the Hamiltonian 
of this theory does not require complex nonlinear terms that characterize many other 
theories or models of glass transitions. 



3. Mean field toy model 

3.1. Model 

Our toy model is a set of oscillators with linear and random nonlinear couplings 
expressed by the following Langevin equations for the iV-component density-like 
variables ai{t) with i = 1,2,- • •, and the M-component velocity- like variables 
with a = z, 2, • • • , M . Here and after we will use Roman indices for the components 
of a and Greek for those of h: 

hi — Kj^oJjo^ -\- Jij(^o,jb{^ (-^^) 

ba = -iba - Ld'^Kjattj '^Jijai^^'^liaj - TSi-j) + fa (13) 

V N 

< fait) >= 0, < fa{t)fp{t') 2-fTSapS{t - t') (14) 

where the summation is implied for repeated indices and overdots denote time 
derivatives. Here 7 is the decay rate of the velocity-like variables 6q and to gives a 
measure of the frequencies of oscillations of the density- like variables Uj . The thermal 
noises fa (t) are independent Gaussian random variables with zero mean and variance 
2jT^ T being the temperature of the heat bath with which the system has a thermal 
contact. The choice of this variance guarantees the proper equilibration of the variables 
{b}. The N xM matrix Kia plays an important role in the model and for later purpose 
we impose the (one-sided) orthogonality 

KiaKi/S = Sap, KiaKja ^ 5i] (15) 

where the last equation is due to the inequality M < N. For M ~ N we can impose 
an additional condition Kia — Sia and hence trivially KiaKja = 5ij. We also note 
that Kia governs linearized reversible dynamics of the model with the dynamical 
matrix ft given by flij = uP'KiaKja- The mode coupling coefficients Jija are chosen 
to be quenched (time-independent) Gaussian random variables with the following 
properties: 



Jija — 0, 



JijaJkip 



{SikSjl + Sil5jk)SaP + Kii3{KkaSjl + KlaSjk) 
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(16) 



where denotes average over the J's. In constructing this model, we were motivated 
by the works ^ in which random couphng models involving an infinite component 
order parameter have been shown to be exactly analyzed by mean-field-type concepts. 

Equation ([l^ ) is analogous to the equation of continuity of fluid and ( p^ ) is like 
the equation of motion where the right hand side is like the force acting on a fluid 
element which corresponds to (|l]) and (^) , respectively. We will eventually take N and 
M infinite with the ratio S* = M/N kept finite. 

One can derive from the Langevin equations (l)-(3) the corresponding Fokker- 
Planck equation for the probability distribution function D{{a}, {b}, t) for our variable 
set denoted as {a}, {6} as follows 

dtD{{a},{h},t) = LD{{a},{b},t) (17) 

where the Fokker-Planck operator is given hy L — Lq + Li + Lmc 
with 

1/9 d ^ 

Lmc = —T^Jija ~^ — Luajba + ■^--uj{uj'^aiaj - T5i_ 



It is then easy to show that the equilibrium stationary distribution (i.e., LD^{a, b) = 0) 
is given by 

De{{a},{b}) = cst.e ^,=1^^"^ Z^»=i st",. (j^g) 
where est. is the normalization factor. 

3.2. Analysis and discussion 

We aim at finding the set of five equilibrium time correlation functions defined by 
Ca{t -t') = ^< aj{t)a,{t') >, Cab{t - t') = -^K,^ < aj{t)b^{t') >, 



j^Kj^ < b^{t)aj{t') >, Cb(t - t') = i 
1 



cut - t') = —K,^ < b^{t)a,(t') >, Cb{t ~t') = — < b^{t)b^{t') >, (20) 



Cf (i - t') = —K.^Kj^ < a,{t)a^{t') > 

It turns out that we need to have the last correlation function to close the self- 
consistent set of equations for the correlators when M < N. Note that for the case 
M = N, we can take = Sia, and then C^ {t - t') = Ca{t - t'). 

In order to obtain this self-consistent set of equations, it is most convenient to 
adapt the generating functional method from which one can write down the set of 



effective linear Langevin equations valid in the limit of Af, N ^ 00. We refer |17 for 
further details. ^From this effective Langevin equations, one can readily derive the 
following closed self-consistent equations for t > for the five correlators: 

Cait) = d*Cba{t) - ^aa ® Ca{t) - 5*^ab ® Cba{t) (21) 

Cba{t) = ~lCba{t) - w'Cf (0 - Eba ® Cf (0 - Ebb ® Cba{t) (22) 

Cabit) = Cbit) - ^aa ® Cab{t) " ^ab ® Cb{t), (23) 

Cb{t) = -lCb{t) - iO^Cabit) - ^ba ® Cabit) - ^bb <E) Cb{t) (24) 

C^{t)=Cba{t)-J:aa^C^{t)-J:ab^Cba{t) (25) 
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where, for any function X{t), X ® a{t) = dt' X{t - t')a{t'). The equations (|2^)- 
( p5| ) constitute the self-consistent equations for the 5 correlators Ca{t), Cba{t), Cab{t), 
Cb{t), and C^{t). This set of equations can be solved numerically with the initial 
conditions Ca(0) = (0) = T/w^ Ca6(0) = C6a(0) = 0, and Cb(0) = T. Here the 
kernels E's are given by 

T.aa{t - t') = 5*^{Ca{t - t')Cb{t - t') + 5*Cab{t - t')Cba{t - t')) , 
E„fc(t - t') = -2S*^Cait - t')Cba{t - t') 
Y.ba{t - t') = ^28*-^Ca{t ~ t')Cab{t - t' ) , 

These kernels arise from the non-linear mode coupling terms in and (^8|). Note 
that the correlator C^{t, t') is not involved in the S's. 

For further analyses it is very convenient to work with the equations of the Laplace 
transformed correlation functions defined as C^{z) = dte'^*^ C{t). Performing 
the Laplace transformation of the self-consistent equations and rearranging them we 
obtain 

T 1 



(27) 



r<L ^ ^(1 ~ ^afc(^)) /28^ 



' (z + ELW)(z + 7 + Sb6W)+c.2(l-S^,W)2 ^ ' 

T{z + S"^ (z)) 
(z + + 7 + S Jz)) + - 



^^(-)- .,.....,....,. "^n"^!:^.. ...... (30) 



(31) 



z + 7 + i;6(,(z) 

For (5* = 1 where M = N and iiTia = ^i^, C^{z) = C^^{z) reproduces the 
equation derived in , apart from the wave number dependence. Note that if we 
put S^^(z) = S^(j(z) — by hand, ( |2^ ) or ( |3l| ) gives a closed equation for Ca(i) 
alone. This equation is nothing but the Leutheusser's schematic MC equation giving 
a dynamic transition from ergodic phase to nonergodic one [Q. But in reality 
and Eab can not be ignored a priori and our numerical solution strongly indicates 
that the system remains ergodic for all temperatures due to the strong contribution of 
these so called hopping terms. Furthermore these hopping terms cannot be made sclf- 
consistently small as temperature is lowered. Therefore the density correlator does not 
show a continuous slowing down with lowering temperature. This result was striking 
to us since usually a mean-field-type theory, such as the dynamics of the spherical 
p-spin model in the limit of ^ oo, often gives a sharp dynamic transition JlGf . 

Thus it is very difficult to understand the idealized MCT without relying upon 
uncontrolled approximations. It is also interesting to note that the ergodicity restoring 
process in our toy model (represented by the kernels T,aa and Eat) has nothing to do 
with a thermally activated energy barrier crossing since the quadratic Hamiltonian in 
our model does not possess such a barrier. 
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4. Reduced Fokker-Planck equation for the density-like variables 

Possibility of nonergodic states in our model can be seen more directly by adiabatically 
eliminating the variables {6} in the limit of large 7 and obtaining the reduced 
Fokker-Planck equation for the distribution function D{{a\,t) containing only the 
{a} variables: 

dD{{a},t) ~ 

— = CFpD[{a},t) 

^ f d 



where we have defined the Fokker-Planck operator Cpp through the second member 
of the above equation. Here the diffusion matrix Qij{{a}) is given by 

Q,j{{a}) = -M,c.Mjc. (33) 
7 

M^a = + —=Jikaak (34) 



The positive semi-definiteness of the diffusion tensor Q is very easy to show because 
for an arbitrary N component real vector we have 

T 2 

XiQtjXj = — ^ {MiaXi) > (35) 

^ a 

The crucial point is that the diffusion matrix Qij is singular for M < N, i.e., 
det|Q| = giving rise to zero eigenvalues for Q [0. This implies that the Fokker- 
Planck equation ( ^ ) can have nonequilibrium stationary solution other than the 
equilibrium one, Z?e({a}) = est. exp(— w^a|/2T). This nonequilibrium stationary 
solutions are precisely the kind of nonergodic states found numerically in the present 
toy model. In fact a class of stationary solutions is given by 

Ds{{a})^J^{^,a,)e-^^^ (36) 

where is an eigenvector of the diffusion matrix Q^j with zero eigenvalue and J-'{x) a 
non-negative function. If the function J-{x) is a constant, then Ds{{a}) ~ De{{a}) is 
the equilibrium distribution, otherwise it is a nonequilibrium stationary distribution. 
Hence the model is nonergodic for < 5* < 1. 

Let us now investigate this point somewhat further. We first define the following 
non-negative ratio: 

R(^{a},t)^SMA (37) 
De{{a}) 

and then introduce a Boltzmann's H-like quantity as 

n{t)= [ d{a}D{{a},t)\nR{{a},t) (38) 



where the integration is over all the variables in the set {a}. Using normalization 
property of the distribution function and integrating by parts assuming a natural 
boundary condition we find 

?i(i) = J d{a}D{{a},t)C''pp{{a})\nR{{a},t) (39) 
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^-(W)-(^-^«.)q.(W)^ (40) 

Now we readily verify that 
Substituting this into ( |39| ) we finaUy find 

d\nR{{a},t) dlnR{{a},t) 
where the contribution of the first term in ( ^ ) vanishes since 



4p({a})lni?({a},t) = -4pi?-Q,,- — -— (41) 



n{t) - - d{a}D{{a], t)Q,,{{a]) (42) 



I 



diaj^C^p ■ ■ ■ = J d{a}D,Clp ■ ■ ■ = j d{a} {CppDe) • • • - (43) 



This result shows a kind of Boltzmann's H-theorem jTsj (See for a discussion for 
general master equations with detailed balance.): 

nit) < (44) 

Let us suppose that we have performed a transformation of the variables {a} 
{s}, D{{a})d{a} D{{s})d{s}, i?({a}) i?({s}) etc. so that the matrix Qiji{a}) 
is diagonalized: 

QviM) ^ Hi^mj (45) 

with the eignevalues Ai({s}) non-negative functions of {s}. For the case of diagonalized 
Qij, (||), the resuhs (|4|) and (|3|) become 

/ d{s} 5: A. i{s})Di{s}, t) (^^^^^^) ' > (46) 

The stationarity condition ^Z)({s}, t) — Q or 'H{t) — then implies 

/ d{s} E ms})Di{s}, t) (^^^^^^) ' - (47) 

If we denote those subset of the variables of {s} with positive eigenvalues as Sq,, a — 
1, 2, • • • , M'(< M), which we denote as {s}', we must have in the region l)({s}, t) > 0, 

^^4^=0f-A„(M)>0 (48) 

We note that validity of the condition (|4^) in general depends on the regions in the 
space of the variables {s} through {s} dependence of the A's. Thus, in a particular 
region in which Aq({s})D({s}) are positive definite, i?({s},t) does not contain s^. If 
we denote the remaining set of the variables {s} than those corresponding to positive 
eigenvalues as u^, /? — M' + 1, M' + 2, ■ ■ ■ , N , which are corresponding to the zero 
eigenvalues, we should have R{{s},t) = R{{s}' , {u} , t) — R{{u},t). Consequently, 
from (^T\), the general form of non-equilibrium stationary state distribution function 
is 

Ds{{s}', {u}) = RsiMmuy, {u}) (49) 
The previous result, (|36|), is a special case of this general form. 
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5. Summary and discussion 

Our toy model analyses show that the so-called hopping processes introduced as a 
correction to the idealized MCT Q are the effects of the velocity-like variables, which 
still enter in the framework of the meanfield treatment, and are not directly related to 
the barrier-crossing mechanism. In our toy model, strength of the hopping processes 
can be adjusted through the parameter 5* . Possible nonequilibrium stationary states 
are connected to the singularity of the diffusion matrix of reduced Fokker-Planck 
equation of Section^ involving only the varaibles {a}. This feature, suitably extended 
to general master equation, would be common to many kinetically constrained glass 
transition models [p^ and is worth further exploration. 
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